Spin polarization of strongly interacting 2D electrons: the role of disorder. 
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In high-mobility silicon MOSFET's, the g*m* inferred indirectly from magnetoconductance and 
magnetoresistance measurements with the assumption that g* hbHs = 2Ef are in surprisingly good 
agreement with g*m* obtained by direct measurement of Shubnikov-de Haas oscillations. The 
enhanced susceptibility x* ^ (<7* m *) exhibits critical behavior of the form \* oc [n — no)~ a . We 
examine the significance of the field scale H s derived from transport measurements, and show that 
this field signals the onset of full spin polarization only in the absence of disorder. Our results 
suggest that disorder becomes increasingly important as the electron density is reduced toward the 
transition. 
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Two dimensional systems of electrons |l]-|3| and holes 
have been the focus of a great deal of attention 
during the last few years. In contrast with expectations 
for noninteracting |Q or weakly interacting ||] electrons 
in two dimensions, these strongly interacting systems 
exhibit metallic behavior in the absence of a magnetic 
field: above some characteristic electron (hole) density, 
n c , their resistivities decrease with decreasing tempera- 
ture. Whether there is a genuine metallic phase and a 
true metal-insulator transition in these materials contin- 
ues to be the subject of lively debate |J. 

Experimental results have been obtained in the 2D sys- 
tem of electrons in silicon MOSFET's that indicate that 
the response to a magnetic field applied in the plane of 
the electrons increases dramatically as the electron den- 
sity is decreased toward n c . Based on a study of the 
scaled magnetoconductance as a function of temperature 
and electron density, Vitkalov et al. jl(|] have identified 
an energy scale A that decreases with decreasing density 
and extrapolates to zero in the limit T — > at a density 
no m the vicinity of n c ; this was interpreted as evidence of 
a quantum phase transition at no. From studies at very 
low temperatures of the magnetoresistance as a function 
of electron density, Shashkin et al. 1 1 inferred that the 
two-dimensional system of electrons in silicon inversion 
layers approaches a ferromagnetic instability at the criti- 
cal density n c for the zero-field metal-insulator transition. 
From a determination of the enhanced spin susceptibil- 
ity derived from Shubnikov-de Haas measurements down 
to low densities, Pudalov et al. [Q have claimed there 
is no spontaneous spin polarization for electron densities 
above n = 8.34 x 10 10 cm~ 2 « n c , although they could 
not exclude this for lower densities. The possibility that 
a magnetically ordered phase exists in the limit T — ► in 
dilute two-dimensional silicon inversion layers is intrigu- 
ing and bears further investigation. 

In this paper we show that there is very good agree- 



ment between values reported for g*m* as a function 
of electron density in high-mobility silicon MOSFET's 
obtained directly from measurements of the Shubnikov- 
de Haas oscillations [jl2) and those inferred indirectly 
from magnetoconductance and magnetoresistance mea- 
surements by two different groups using different meth- 
ods of analysis and the assumption that g* [XbH s = 2Ep 
p0| , pd 13 1 . Here g* is the enhanced g-factor, m* is the 
enhanced electron mass, us is Boltzmann's factor, Ep is 
the Fermi energy, and H s is a characteristic field scale 
determined by different methods from in-plane magne- 
toconductance jl(| and magnetoresistance |ll|] experi- 
ments. The enhanced susceptibility x* 0(1 (s* 771 *) ex ~ 
hibits critical behavior of the form \* °^ (n— no)~ Q . Data 
from the three experimental groups yield exponents a of 
0.23, 0.24 and 0.27, and critical densities between 0.88 
and 1.04 x 10 11 cm~ 2 . We examine the significance of the 
field scale H s , and show that this field signals the onset 
of full spin polarization only in the absence of disorder. 
Our results suggest that disorder becomes increasingly 
important as the electron density is reduced toward the 
transition. 

Measurements of Shubnikov-de Haas oscillations in 
high-mobility silicon MOSFET's with high electron den- 
sities have shown that the magnetic field required to 
achieve complete polarization of the electron spins is ap- 
proximately the same as that required to saturate the 
magnetoresistance to a constant value 1 16 - 18 1 . For the 
relatively high densities used in these experiments, the 
field H p correponding to saturation of the magnetoresis- 
tance is approximately the same as the field H a above 
which there is apparent saturation of the magnetocon- 
ductance. As we show below, this equivalence breaks 
down at lower densities. A clear example is illustrated in 
Fig. [I], where the resistivity and conductivity are shown 
as a function of in-plane magnetic field for a silicon MOS- 
FET with electron density near the critical density, n c , 



1 



for the metal-insulator transition. The saturation field 
Hp derived from the resistivity is considerably larger than 
the field H a above which the conductivity saturates. This 
can be understood with reference to the band diagrams 
shown as insets to Fig. [I]. In the absence of disorder, 




FIG. 1. For a silicon MOSFET with electron 
density 0.94 x 10 11 cm~ 2 , the conductivity (left 
curve) and resistivity (right curve) are shown as 
a function of in-plane magnetic field at temper- 
ature T = 0.26 K; the saturation fields deduced 
from the resistivity and the conductivity are la- 
beled Hp and H ai respectively. The insets show 
schematic diagrams of the electron bands (see 
text for discussion). 



all electron states are extended, band-tailing plays a neg- 
ligible role, and full spin polarization is achieved when the 
Zeeman energy is sufficient to completely depopulate the 
minority spin band: 



g fisHband — 2E F : 



(1) 



here g* is the enhanced g-factor, /is is Boltzmann's fac- 
tor, Hband is the magnetic field required to fully polarize 
the system in the absence of disorder, and Ep is the 
Fermi energy. Disorder is weak at high electron densities 



and one expects H t 



band 



Hn ~ H„ 



tails become more important, complete spin alignment 
requires the application of a larger magnetic field to fully 
polarize the tail states as well as the extended states: 



g l^BHtail+ba 



2Ef 



(2) 



where we've assumed the band tail has an effective energy 
width S JT|. 

Except very near the transition, the number of states in 
the band tails in the case of samples of reasonably high 
mobility is much smaller than the number of extended 
states; at the same time, the energy width S becomes 
appreciable as the density decreases and the disorder in- 
creases. The field required to align the electrons in the 
higher mobility band states can thus differ substantially 
from the magnetic field needed to polarize all the elec- 
trons |H| . While the (small number) of tail states make 
a minor contribution to the conductivity, the resistivity 
is considerably more sensitive to the low-mobility states 
in the tail of the distribution, and consequently H p > H a 
as is evident in Fig]l]. We suggest that H a w H oan d and 
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FIG. 2. The fractional difference (H p — H a )/H a (open 
symbols) and 2/er (closed symbols) versus electron den- 
sity; Hp and H a are the saturation fields deduced from 
resistivity and conductivity curves, respectively. The 
open circles and closed symbols refer to data taken at 
0.26 K. The open squares are data obtained at 0.1 K 
on a different MOSFET. 



The fractional difference between H p and H a , 
AH/H = (Hp — H a )jH a , is shown as a function of elec- 
tron density in Fig. ^; AH/H increases rapidly with 
decreasing electron density when disorder becomes more 
dominant. The quantity 2/ a is plotted for compari- 
son through the following argument. For weak scat- 
tering, the parameter S is on the order of the scatter- 
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ing rate: 6 ~ Ti/t. With Eqs. [I] and ||, this gives 
AH/H = 5/2E F = h/2{E F r). Using the expression for 
the Drude conductivity a = ne 2 T/m* 1 and the Fermi 
energy Ef/H — (nh)/ g v g s m* with a valley degener- 
acy g v = 2 and spin degeneracy g s = 2, one obtains 
AH/H = (e 2 /h)(2/a). The correlation between AH/H 
and 2/ a is evident in Fig. ||. 

In an earlier paper [fTo|| , we showed that the magneto- 
conductance of silicon MOSFET's can be scaled onto a 
single curve by plotting [a(H)—a(Q)]/[a(H — oo) — c(0)] 
as a function of H/H s . The parameter H s obtained by 
this method is proportional to H a discussed above. For 
high densities where disorder plays a small role, the mag- 
netic field H a needed to saturate the conductivity is very 
nearly equal to the field required to obtain full spin po- 
larization. At lower densities, the saturation fields de- 
duced from the resistivity and the conductivity are not 
the same, and we have argued that the difference is as- 
sociated with the effect of electrons in the states in the 
band tails. We've suggested that H a is the magnetic field 
required to polarize the band states; the Zeeman energy 
and g*m* are then given by Eq. [I] with H band = H a . 
The tail states remain unpolarized in H = H a . However, 
except perhaps very near the transition (or in samples of 
very low mobility), they represent a small fraction of the 
electrons, so that the system is close to full spin polar- 
ization. 
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FIG. 3. The inverse of the enhanced susceptibility 
Xo/x* versus electron density obtained by Vitkalov 
et al. [nj, Shashkin et al. H], and Pudalov et al. 
P2| . Data are normalized to the Shubnikov-de Haas 
values at high densities. The curve is a fit to the 
critical form xa/x* = A(n — n^) a for the data of 
ref. |1(J (excluding the point shown at Xo/x* = 0)- 



Fig. |^ shows 2mo/m*g* — Xa/x* as a function of elec- 
tron density n s obtained from our data |f0[ , by Sashkin 
et al. Q, and Pudalov et al. Q. Here x*/Xo is 
the enhanced susceptibility normalized to its free elec- 
tron value, and Xo/x* is its inverse. The closed cir- 
cles denote values obtained from scaling our data for 
the in-plane magnetoconductance and the assumption 
that g*HBH a = 2Ep\ the open circles were obtained 



by Shashkin et al. If |from magnetoresistance measure- 
ments using a different data-fitting procedure and the 
same assumption as above; the squares are from direct 
Shubnikov-de Haas measurements of Pudalov et al. [[l2) . 
The data of Shashkin et al. decrease somewhat more 
rapidly at low densities than the others. However, the 
three sets obtained by different groups using different 
measurements and different methods of analysis agree 
surprisingly well. Again, this indicates that the small 
number of states in the band tails in high-mobility MOS- 
FET's play a neglibible role. A fit to the critical form 



Xo/x* oc(n- n ) a , 



(3) 



2 ; for the mag- 
pll a = 0.27, 



yields the following values for the three data sets con- 
sidered: for the Shubnikov-de Haas data of Pudalov et 
al. HI a = 0.23, n = 0.96 x I0 11 cm" 
netoconductance data of Shashkin et al 
n = 1.04 x I0 11 cm~ 2 ; and for our data |H| a = 0.24, 
n = 0.88 x I0 11 cm" 2 . 

We have argued above that for high- mobility samples, 
the difference (H p — H a ) is associated with the effect of 
a small fraction of the electrons in the band tails. The 
characteristic field H s obtained in our earlier work was 
determined from scaling the magnetoconductance, which 
is a measure of the field required to align the band states 
while leaving a few electrons in the tail states unpolar- 
ized. Shashkin et al. determined a field scale by match- 
ing magnetoresistance data at low magnetic fields; close 
examination shows that this procedure does not produce 
a match at high fields (note that their data is shown on 
a logarithmic scale, which deemphasizes differences be- 
tween the curves at high values of magnetic field) . Both 
methods are sensitive to the contribution of the extended 
state and minimize the effect of the states in the band 
tails. These procedures yield reliable measures for the 
behavior of the system at high electron densities where 
disorder does not play an important role. This accounts 
for the surprisingly good agreement between the g*m* 
obtained from transport experiments and those found 
by direct measurement of the Shubnikov-de Haas oscilla- 
tions. At densities very near the transition (and for very 
low mobility MOSFET's) one should expect this corre- 
spondence to break down as disorder becomes more dom- 
inant. We suggest that an understanding of any phase 
transition that occurs in this regime must incorporate the 
effect of disorder in a central way. 
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